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PVagmentation Production of Triply Heavy Baryons at the LHC 

M.A. Gomshi NobarjO and R. Sepahvancjj 
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The triply heavy baryons in the standard model formed in direct c and b quark fragmentation are 
the Q.CCC, ^ccb, 0.cbb and Q,bbb baryons. We calculate their fragmentation functions in leading order of 
perturbative QCD. The universal fragmentation probabilities fall within the range of 10~^ — lO'^^.We 
also evaluate their cross section at the LHC {^/s = 14 TeV) using next-to-leading order matrix 
elements for heavy quark-antiquark pair production. We present the differential cross sections as 
functions of the transverse momentum as well as the total cross sections. They range from a few nb 
to a few pb. 

PACS numbers: 13.87.Fh, 13.85.Ni, 12.39.Hg 



I. INTRODUCTION 

Heavy hadrons have been the focus of attention due to 
their interesting properties. The study of production and 
decay of such particles is interesting in two aspects. In 
the first place the question is whether QCD is the right 
theory to predict the properties of such objects through 
confirmation of the standard model predictions with ex- 
perimental data. Secondly this research investigates the 
basic properties of the weak interactions at the funda- 
mental level. These states have in general a large number 
of decay modes so that their observation and measure- 
ment of their properties require a large number of them 
to be produced. Their cross section at e+e" collisions is 
very small therefore their identification needs a messier 
environment of the hadronic collider. 

In the framework of the quark model, heavy baryons 
fall into three categories. States containing one heavy fla- 
vor such as Ac and Af,, are interesting states due to the 
fact that they carry the original heavy flavor polariza- 
tion [1]. Their production has been studied in interest- 
ing models [2]. They are also being studied experimen- 
tally [3] . The second category involves baryons with two 
heavy flavor like the states 'E.cc, "^bb and Sf,c [4]. They are 
treated within the approximate quark-diquark model [5] . 
The model treats the production of the so called diquark 
perturbatively similar to B^ states [6]. Then, it can be 
proved that the formation of a baryon out of a diquark 
is almost the same as the fragmentation of an antiquark 
into a meson [7]. In this way one obtains the fragmen- 
tation functions, the total production probabilities and 
their event rates in a desired collider. Indeed the light 
degree of freedom within these states does not allow full 
perturbative calculation. 

In the third category, we have baryons with three heavy 
constituents. Since the top quark cannot take part in 
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strong interactions [8] , there remains only the charm and 
bottom quarks to form such baryons. There have been 
attempts to evaluate the production of ^ccc and ^bbb iu 
e''"e~ and hadron colliders in the quark-diquark model [9] 
and also using perturbative QCD [10]. The results from 
e'^e~ annihilation are very small indeed [11]. However 
sizable rates are expected in energetic hadron colliders 
[12]. Therefore the standard model production rates of 
these bound states can be compared with experimental 
data [13]. 

Consistent with the quark model of hadrons, the spec- 
troscopy and production mechanism of heavy meson and 
baryon states have been treated satisfactorily. Specially 
the hadrons which contain c and b quarks or anti-quarks, 
are accounted for in the heavy quark limit where the 
hadronic bound state is understood and the perturbation 
theory is applied for the process of their production. This 
has been successful in the treatment of B^ states both in 
theory [6], and in experiment [14] and also in the produc- 
tion of heavy diquarks in the treatment of doubly heavy 
baryons [4] . In this work we shall apply this procedure to 
the case of triply heavy baryons to obtain their fragmen- 
tation properties and cross section at the LHC. Many of 
these states may be observed at existing hadron colliders, 
specially at the Tevatron, however some others have very 
low event rates. Therefore we have chosen the LHC for 
the sake of integrity. Therefore in this work we consider 
a framework which treats all triply heavy baryons and 
obtain their fragmentation functions and estimate their 
production at the LHC. 

Our plan is as follows. In section II we provide a gen- 
eral discussion of the fragmentation process of S-wave 
triply heavy baryons and calculation of their fragmenta- 
tion functions. In sections III and IV we calculate the 
fragmentation functions for c — > i^cch and h -^ ^ccb which 
we have chosen to be the basic ones such that the other 
functions could be obtained from them by appropriate 
choices of quark masses and other baryon characteristic 
parameters. The inclusive production of these states at 
the LHC is studied in section V. Finally we discuss our 
results in section VI. 
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FIG. 1: The lowest order Feynman diagrams contributing to 
the fragmentation of a heavy quark (Q) into a triply heavy 
baryon (B). The four momenta are labelled. 



II. FRAGMENTATION OF TRIPLY HEAVY 
BARYONS 

The fragmentation functions are process independent 
and can be applied to the e+e~, partonic and hadronic 
production processes. At sufficiently large transverse mo- 
menta, the dominant production mechanism is actually 
the fragmentation, the production of a parton with high 
transverse momentum which subsequently splits into a 
triply heavy state and other partons. Fig.l shows the 
fragmentation of a heavy quark Q into a triply heavy 
baryon B{QQ'Q") in lowest order perturbation theory. 
We will calculate such Feynman diagrams. 

The fragmentation of a parton into a baryon state is 
described by fragmentation function D{z,fx), where z is 
the longitudinal momentum fraction of the baryon state 
and /^ is a fragmentation scale. The fragmentation func- 
tion for the production of an S-wave triply heavy baryon 
B in the fragmentation of a quark Q is obtained from 
[15] 



X<53(p + s' + t'-p)d3pd3s'd3t', (1) 
where four momenta are as labelled in Fig 1. Tb is the 



amplitude of the baryon production which involves the 
hard scattering amplitude Th and the non-perturbative 
smearing of the bound state. The average over initial and 
the sum over final spin states are assumed. 

To absorb the soft behavior of the bound state into 
the hard scattering amplitude we have used the scheme 
introduced in [16]. The probability amplitude at large 
momentum transfer factorizes into a convolution of the 
hard-scattering amplitude Tu and baryon distribution 
amplitude 0b, i.e. 



Tb(p,p,s' ,t')^ [dx]TH{p,p,s',t',Xi)(j)B{xi,w'^). (2) 

where 4'b is the probability amplitude to find quarks co- 
linear up to a scale of w'^ in the baryonic bound state. 
In (2), \dx] — dxidx2dx2S{l — xi — a:;2 — x^) and Xi's 
are the momentum fractions carried by the constituent 
quarks and finally Th is written in the following form in 
the old fashioned perturbation theory to keep the initial 
heavy quark on mass shell. 



Th = 



mir^as {2mQ, )a, {2mQ„ )MCf 



2^2poP,s'X 



(Po + S'o+t'o -Po) 



.(3) 



Here F represents an appropriate combination of the 
propagators and the spinorial parts of the amplitude. 
as = <7^/47r is the strong interaction coupling constant. 
Cf is the color factor and M = mmQmQ>mQ>/ with m 
being the baryon mass. Since we ignore the virtual mo- 
tion of the baryon constituents, we propose a delta func- 
tion type function to represent the probability amplitude 
of the baryon state 



4>B = Ib S< Xi 



rrii 



(4) 



where fs is the baryon decay constant and is introduced 
similar to meson decay constant Jm- Putting this ex- 
pression and (3) in (2) and carrying out the necessary 
integrations, wc find 



Tr 



l&TT^as{2mQ,)as{2mQn)MCF.fE 
2^2po%s'X 
F 



{Po +s'o +t'o-Po)' 



(5) 



With this amplitude wc find the fragmentation function 
as 



D^{z,fi) ^ 32[Tr^a,{2mQ,)a,{2mQ„)MCFfE 



iErr5'(p + s' + t' 



PoPoS'ot'oiPo + S'o + K - Po} 



'^' ;A^pd?s'dH'. 



r 



(6) 



To proceed we need to specify our kinematics. 



We let the baryon move in the z direction after produc- 
tion, neglecting the virtual motion of the constituents. 
The initial state heavy quark has a transverse momen- 
tum which should be carried by the two antiquarks away 
through the final state jet. We have assumed that there 
will be only one jet in the final state. This assumption 
is justified due to the fact that the very high momentum 
of the initial heavy quark will predominantly be carried 
in the forward direction. Due to momentum conserva- 
tion, the total transverse momentum of the two jets will 
be identical to the transverse momentum of the initial 
heavy quark. Therefore the antiquark's contributions to 
this jet are assumed to be proportional to their masses. 

The fragmentation parameter z, is defined as usual, i.e. 



{E + P\\)b 
{E+p«)q-' 



(7) 



which reduces to the following in the infinite momentum 
frame which we have adopted for our study 



Eb 



(8) 



Now we set up our kinematics. According to Fig. 1 
the baryon takes a fraction z of the initial heavy quark's 
energy (each constituent a fraction of xi, X2 and x^) and 
the two anti-quarks take the remaining I ~ z (x4 and 
X5 each). Thus the four momenta of the particles are 
parameterized as 



p^ — Zpo So — XlZpo To = X2Zpo 
s'o = X^{l-z)po i'o == 2:5(1 - z)po 



to 



X3ZP0 



(9) 



where the condition xi + X2 + xt, = 
regarding our assumptions we have 



1 holds. Moreover 



s T = x^pt 



t'r 



X5PT, 



(10) 



along with the constraint of 2:4 + 0:5 = 1 . 

Fig. 2 shows the lowest order Feynman diagrams for 
the fragmentation of flccc and Qbhb (a-jb), i~lccb in the c 
quark fragmentation (c,d) and Qccb in the b quark frag- 
mentation (e,f). There are similar diagrams contributing 
to flcbb fragmentation in b and c quark fragmentation 
which are simply obtained by interchanging the c and b 
quarks in (c,d) and (e,f ) respectively. Let us first consider 
the case of c ^- flccb- 
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FIG. 2: The lowest order Feynman diagrams which con- 
tribute to different Q baryon production. While (a) and (b) 
shows c and b quark fragmentation into ^ccc and Qtbt re- 
spectively, (c) and (d)contribute to fleet production in c and 
(e) and (f) in b quark fragmentation to the same state. In- 
terchange of c < > b in last two pairs give the contributing 

diagrams to Qcbt in 6 and flcbb in c quark fragmentation re- 
spectively. 



III. Qccb IN c QUARK FRAGMENTATION 

Here the diagrams (c) and (d) in Fig 2 are relevant. In 
each diagram there are three propagators. In this specific 
case for the diagram (c) we find the combination 



G 



1 



8to^TO2to'/'^(z) ' 



(11) 



where one third of the power of / comes from each prop- 
agator, m' — (jTii + 7712) and that 



m 1 ~ z to' / Pt \ ^ 
■'^'"'^ ^ 2to' z ^2^V-^'^){1-z)' 



In the case of diagram (d) we have 
G' 



(12) 



(13) 



Idrnfm^f^iz) 
We put the dot products of the relevant four vectors 



in the following form 

I 2 1 a I 2 

s .r — ra^a r.p — rriip s .p ^ rriij 

where 

1 



m 1 — z m' / P% \ z 

— — + — 1 + ^ r- 

m z m \ m'-' /I — 



'-^ 



m \ mi / mi z 



(14) 



(15) 



h(i+^l)(i- 

m' \ mf/ 
^' (.,Pt\ 1 


-^) 


mi \ 


m'V 1- 


z 



2p: 



mim 



(17) 



In obtaining the above results we have used (9) and (10) 
with xi = a;2 = mi/m;x3 = m,2/m,, Xi = m,i/m,' and 
Xs — m,2/m,'. 

(16) In this case the F in (6) reads 



r= Guit)Yvit') u{r)-fp{^ + mi)-f^u{p) u{s)-t''v{s') + G'uis)-t''vis') uir)-f^i4' + mi)-fpu{p) u{t)Yv{t'). (18) 



From which we find 



where 



FF = 



G^ [(/ - mi)7^(^ + mi)Y\ T^^p. [{f - m2)7"(|^+ ^2)7" 
G'2 [if - m2)7"(|^+ ^2)7"] T,^,p [(/ - mi)7^(^ + 7^)7'^ 
GG' [(/ - mi)7^(^ + 7711)7"^] T^,,p [if - TO2)7"(/+ "^2)7" 
GG' \{f - m2)i%i+ ^2)7"! r.^p. f(/ - mi)7^(^ + 7711)7'^ , 



Tpapu = {i> + mi)-fp{i + miYf^lf + miYfpli + mi)^^, 



(19) 



(20) 



Tap^^p = i^ + mi)-fcr{i' + mi)jp{f + mi)7i.(/ + 7771)7 



P' 



(21) 



TcTf^piy ^ {i> + ™i)7m(/ + ^lYfaif + mi)j^{^ + mi)7 



P' 



(22) 



Tcri^piy ^ {i> + mi)-f„{i + mi)j^{f + 7771 )7p(/ + mi)j^ 

I 



(23) 



Next we consider the phase space integrations. Note that and 



/ = 



where 



6^p + s'- 



t'-p) 



d^p = 



Po 



Po{Po +ko+K- Pof [7777r7'5(z)] 



r, (24) 



p% ml 

mm m z 
Pt 



m 



('-!&) 



1 



'2/1 • (25) 

777 \ 777^ / 1 — Z 

Here instead of performing transverse momentum inte- 
grations we replace the integration variable by its average 
value in each case. Therefore we write 



/ 771/^ „/ ^J3„/ _ / 



TT^/^ „/ \ J„' j2„ 



j F'{z,t'T)dh' = j 



F'{z,t'T)dt'i^d% 
mlt',F'{z, {t'j.')) 

2,1 w r ™2 , 2 

'm2toF z,^{pT 

L 777 



(27) 



Putting all this together back in (6), we obtain the 



frA.O'vnenta.tinn fnnr'tinn for 



0.__t as fnl1n-wfi 



Dr 



,(2, Mo) 






6 - /3 - 7 + ^^(l + 7) + 4a\a - 1) -2 + a(/3 - 1) + /3 + 7 



+^2(4 + 5(3 + 37) + a(7 + 6/3 + 57) + 8a3 5 + /3 + 0^(6 + 9/3) + a{7 + 7/3 + 7) 



-2a 



66 + 23/3 + 277 + a^{l + 7) + a^(52 + 37/3 + I97) + a(115 + 70/3 + 577) 
182 + 71/3 + 7I7 + 5a-''(l + 7) + ^^(168 + 141/3 + 477) + a(311 + 198/3 + I337) 

I 



(28) 



Here f{z) given by (12) is due to the propagators and 
g{z)coines from the energy denominator (25). a, (3 and 
7 are for dot products given by (15)-(17). We have set 
a = mi/m2. 

It is clear that the interchange of c <-^ fe in the above 
function will provide the fragmentation function for b -^ 
ricbb in agreement with our direct calculation. 



IV. noo6 IN b QUARK FRAGMENTATION 

Now let us consider the process of 6 — > ^ccb- In this 
case regarding the diagrams (e) and (f) in Fig. 2 and 
using the above procedure we find for the propagators 



G^G' 



1 



8mfm'/'^(2)' 



where f'{z) is 
/'(^)-l + 



m 1 — z mi 



Am 



1 z 



1Pt\ z 
Am\) I- z' 



(29) 



(30) 



The dot products of the relevant four vectors are put 



in the following form 

/ 21 2 a// 2 / 

s .r — m2a r.p — m2P s .p ~ m^'y 

where 



a 



m 1 — z 
2to2 z 



-I 



mm2 



2z 



Ami) I- 



P' 



TO2 



mi J 



m 1 

7712 Z 



^ ^ 2 



i(i^)Vi+4)(i-z) 

2\mi/ V 7712/ 



4 ?Tl^ / 1 — z 



(31) 



(32) 



(33) 



(34) 



Note that the x's in (9) read as xi ^ x:} = mi/m, X2 — 
TO2/TO and X4 = X5 = 1/2 in this case. Here the T in (6) 
has the following form 



Therefore 



r' ^ Gu{t)Yv{t') u{r)-fp{i + m2)-fMp) uish^vis') + G'u{s)-f''vis') u{r)-f^{i' + m2hpu{p) u{t)Yv{t'). (35) 



where 
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^ fiapiy 
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-7771)7^(^ + 7771)7'^ 
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-777l)7"(^ + 7771)7"^" 
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TL^p ^ i^ + m2hai4' + m2hJf + m2h^{4' + m2hn, 



(36) 

(37) 
(38) 



Kt^pu = {i> + TTi2YlM + 'm-2haif + ?T^2)7l.(^ + n^2)7p, 



TL. 



{i) + m2ha{4 + 'm2)^p{f + m2)"fp{i' + m2)ju- 



(39) 

(40) 



Finally similar to our previous treatment of c ^ ^ccb, we obtain the fragmentation function for b -^ flccb as 



^4a,(2mi)V|C| 



64a6TO'2z/'6(z)g'2(z) 

+14a^7' + 20^^(12 + 5/3' + 7' + 12a'7') + a^ 

^ ^20 + 6/3' - 7' + lla'27' + 4a'(ll + 5/3' + 37') | 



ha'^il + 3/3') + aa' -5 + a'^ - 4/3' + 2a'(8 + 9/3') 

2 - /3' + a'2(21 + 10/3') + a'(36 + 34/3' - 37') 



+a'37' 



(41) 



Here g'{z) comes from the energy denominator which in 
this case reads as 



g'i^) 



Pt 



m 1 



21717711 27171 Z 

2toi /_, p^ 



A7n\J 1 — z 



(42) 



r 



Again in this case the interchange of c < — > h will pro- 
vide the fragmentation function for c -^ i^cbb- This also 
agrees with our direct calculation. 

In the equal mass case where a' = a, (3' = /3, 7' = 7, 
f'{z) = f{z) and g'{z) = g{z) the fragmentation function 
takes the form 



Dq^dQQQiz^^io) 



^^as{27nQYflCl 
256mlz,p{z)g^{z) 

X Ue + 15/3 + 157 + a^{l + 7) + a^(40 + 37/3 + II7 + a(75 + 50/3 + 337)|, (43) 



where Q may be assumed to be a c or 5 quark with ttiq 
being respective quark mass. 

The input for the fragmentation functions (28), (41) 
and (43) are quark masses, baryon decay constants and 
the color factor. We have set toi = TOc=1.25 GeV and 
TO2 = TO&=4.25 GeV. For the decay constant and the 
color factor we have taken fB—0-25 GeV and Cp = 7/6 
for all cases. The later being calculated using color line 
counting rule. We have also taken (p^) — ! GeV which is 
an optimum value for this quantity. 



r 



factorization. Essentially this idea incorporates the short 
distance high energy parton production and the long dis- 
tance fragmentation process. Here it is assumed that 
at high transverse momentum the inclusive production 
of triply heavy baryons is factorized into convolution of 
parton distribution functions, bare cross section of the 
initiating heavy quark and the fragmentation function, 
i.e. 



V. INCLUSIVE PRODUCTION CROSS 
SECTION 

Theoretical calculations of the production cross section 
in high energy hadron collisions are based on the idea of 
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FIG. 3; Born diagrams contributing to the calculation of 
heavy quark pair production cross section. 



The physical production rates calculated at all orders 
in perturbation theory would be independent of normal- 
ization/ factorization scale. Such results are not avail- 
able. So that the production cross sections do depend to 
a certain degree on the choices of /i. We will estimate 
the dependence on /i by choosing the transverse mass of 
the heavy quark as our central choice of scale defined by 



M-R = A/pT^(parton) + mg^ 



(46) 



and vary it appropriately to the fragmentation scale of 
our particles. This choice of scale, which is of the or- 
der of pt (parton), avoids the large logarithms in the 
process of the form ln(mQ//L() or ln(pT//^)- However, we 
have to sum up the logarithms of order of fin/mQ in the 
fragmentation functions. But this can be implemented 
by evolving the fragmentation functions by the Altarelli- 
Parisi equation. This equation reads as 




(a) 




(b) ^ 







dy 



tJ--gj^DQ^Hiz,^i) 



Pq^Qiz/y-, n)DQ^H{.y-, m)- 



(47) 



Here the functions D{z^ fi) at the initial scale /Xo are given 
by (28), (41) and (43). Pq-,q{x = z/y, /x) is the Altarelli- 
Parisi splitting function and at the leading order in Ug 
reads 



FIG. 4: Examples of Feynman diagrams up to order of ai 
contributing to quark-antiquark pair production. 
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A:as{ii) f 1 + X 



37r 



1 



(48) 



mentation of the produced heavy quark into a triply 
heavy baryon. Note that the equation is written in the 
factorization scale /i. Furthermore here the factoriza- 
tion, renormalization and fragmentation scales are set to 
be equal. In other words we have set the scale /i in the 
parton distribution functions, subprocess cross sections 
and the fragmentation functions to be the same. Mean- 
while the running scale used in fragmentation functions 
is set to be maximum (/i, /io), where /io is the prescribed 
initial scale for the fragmentation functions. We have 
employed the parameterization due to Martin-Roberts- 
Stiriling (MRS) [17] for parton distribution functions and 
have included the heavy quark production cross section 
up to the order of a^ [18]. The examples of Feynman 
diagrams for the QQ pair production to the order of a^ 
and af are shown in figure 3 and 4 respectively. 

For the LHC the acceptance cuts of pr'^"* > 10 GeV 
and [yl < 1 are chosen where the rapidity y is defined as 



where the running coupling constant as(/z) is evaluated 
at one loop by evolving from the experimental value 
as{Mz) = 0.1172 [19] given by 



as{^J■) 



as{Mz) 



33 -2n 



l + 87r6oas(M2 



ln(^ 



Mz ' 



487r2 



(49) 



Here Uf is number of flavors below the scale ^. The (+) 

prescription reads f{x)^ = f{^) ^ H^ ^ ^) Jn f{x')dx'. 
We note that only Pq^q splitting function appears in 
(47) . This is because the quark Q is assumed to be heavy 
enough to make other contributions namely Pq^g, Pg^Q 
and Pg^g irrelevant. The boundary condition on the evo- 
lution equation (47) is the initial fragmentation function 
Dq^h{z, fJ.) at some scale fi — fXo where its calculation 
is possible. 



VI. RESULTS AND DISCUSSION 



2 ^[E + pL 



(45) 



In the heavy quark limit we have obtained exact an- 
alytical fragmentation functions for S-wave triply heavy 
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FIG. 5: Fragmentation of a c and h quark into possible triply heavy baryons. Note that they are grouped according to their 
heavy contents. Evolution to desired scales are shown for the LHC. We have used two sets of scales p = /^i?/2, 2/i_R, 4p_R (left) 
and n = fiR, 3fj,R, 6fj,R (right). 
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FIG. 6: Differential cross section d(T/dpT[nb/(GeV)] versus transverse momentum pr for production of different triply heavy 
baryons in c and h fragmentation at the LHC. In each graph the distribution is shown for scales specified. Note that the sets 
of scale used in each column of diagrams is the same (the baryons in the left column contain at least two c and in the right 
column at least two h quark). The kinematical cuts imposed are pt"^"* > 10 GeV andlyl < 1. 
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baryons using leading order pcrturbative QCD. The non- 
petubative part of the bound state is treated by employ- 
ing a delta function type distribution function thus ignor- 
ing the respective motion of the constituents. We have 
obtained the fragmentation functions for c -^ i^ccb and 
b -^ flccb at the scale of /io = m + mQi + vfiQn where m 
is the baryon mass and Q' and Q" are specified in Fig. 
1. The functions for c -^ '^cbh and h -^ '^cbh are ob- 
tained simply by the interchange of c <-> 6 in agreement 
with direct calculations. The fragmentation functions for 
c — > ^ccc and h — > f2f,&f) are obtained by setting the c and 
h quark masses to be equal. 

With our choice of quark masses, i.e. mc—1.25 GeV, 
nil,— 4.25 GeV and fip defined by (46) the behavior of 
fragmentation functions as well as the transverse mo- 
mentum distributions of the differential cross sections are 
well analyzed if we put the triply heavy baryons in two 
groups. Those which contain at least two c and those 
at least two b quarks. Therefore while we need to study 
the c -^ ^ccc, c -^ i^ccb, i> — > ^ccb within a lower set of 
^ = A^fl/2, 2)111 and A^r scales, b — > flhbb, b -^ ilcbb, 
c -^ ^cbb would require higher set of /i = ^ir, S^ir and 
6/i7j. This also provides a means by which the sensitivity 
of the results are tested. When the first (out of the three 
selected) scale fi — ij-r/2 (p = ^r) is less than /io, which 
incidently happens for all of our particles, we choose the 
larger of (^, ^o)- 

The behavior of our fragmentation functions along 
with their evolutions aX ji — jjlr/2 (hr), 2^r (3/i and 
4/i/f (6^) using the Altarelli-Parisi evolution equation 
(47) are shown in figure 5 for different Vt states. Note 
that the scales used here in each diagram are the same 
as the ones which are employed in pT distribution dia- 
grams in Fig. 6. Also note that each column of diagrams 
are sketched in separate set of scales. The universal frag- 
mentation probabilities and the average fragmentation 
parameters at /io are shown in table I. The probabilities 
at this table indicate that while some of the states would 
have considerable event rates at existing colliders, others 
are less probable. Therefore here we present their cross 
sections at the LHC at y^=14 TeV. 

The differential cross sections are shown in figure 6. 
The slow fall off of the distributions is expected in the 
framework of our study. We are dealing with a collider 
with large s/s and production of heaviest hadrons in the 
standard model both against a sharp fall off. A look at 
this figure reveals that firstly the differential cross sec- 
tions are sensitive for different scales chosen only at high 
transverse momentum for nearly all states, secondly as 
the number of the b quark is increased in the state, the 
differential cross section is less sensitive for higher scales. 
It is seen in the cases of 6 ^ ^666, b —^ flcbb and c — > Qcbb- 
This was the main reason to choose two different sets of 
scales here. Sometimes it happens that the distribution 
for two different scales cross each other. This occurs for 
b —^ ^bbb at pt=20 GcV and for b -^ ^cbb at nearly 
Pt=15 GeV and 6 -^ ^ccb at pT—20 GeV all for the first 
two scales. This means that rate of decrease in differen- 



Process 


P.P. 


C ^ n^cc 


2.789 X 10"'' 


c -^ r^ccfc 


2.475 X 10~^ 


b -^ n^cb 


2.183 X 10"^ 


b -^ Qttb 


6.459 X 10-^ 


b -^ i^cbb 


5.290 X 10"** 


C -^ flcbb 


1.086 X 10"^ 



TABLE L The universal fragmentation probabilities (P.P.) 
and the average fragmentation parameter (z) at fragmenta- 
tion scale for different Q states in possible c and b quark frag- 
mentation. 



(^)(m°) 

0.521 
0.490 
0.634 
0.534 
0.562 
0.482 



tial cross sections are different in the specified px region 
for the two scales. 

The total cross sections are listed in table II for the 
chosen scales. They range from a few nb to a few pb. 
The decimal places are not realy significant. They are 
kept only for the matter of comparison. A short look 
at table II reveals that although the total cross section 
for some of the triply heavy baryons are small indeed 
(order of pb) and their production needs energetic hadron 
colliders, some others such as & — > flccb and c -^ ^ccc do 
possess larger cross sections of the order of nb and may 
easily be produced at the Tevatron as well. An interesting 
point in table II is that although the total cross section 
for some of the particles such as c — > flccb and c -^ i^cbb 
increase with increasing /i, but this is not the case for the 
rest. Our investigation shows that this depends on the 
range of /i selected and also on the choice of pT*^"* [20]. 
We have also calculated the ratio a{fl)/a{Q) for different 
cases. The results appear in the last column of table II. 
Our evaluation of charm and bottom cross sections at the 
LHC are 0.25497 mb and 0.46812 mb respectively. 

The fragmentation production of doubly heavy 
baryons studied in [4] by Doncheski et al is interesting 
in relation to our work. First of all the fact that the 
Tevatron gives large cross section for charm production 
is refiected in this work. They have obtained nearly equal 
cross sections for E,cc at the Tevatron and at the LHC. 
However for Ecb and Ebb the cross sections are different. 
They report 430 pb, 215 pb and 16 pb for the Tevatron 
and 470 pb, 490 pb and 36 pb for the LHC respectively 
for these states. Although states which we have studied 
are different, but physically our results are comparable 
with the above. 

We would like at the end discuss the uncertainties of 
our results. The choice of quark masses will not only 
alter the fragmentation probabilities, but also the value 
of /i and values of x at which the parton distribution 
functions are evaluated. This will of course be reflected 
on the total cross sections. We have chosen rric = 1.25 
GeV and mt = 4.25 GeV which are the optimum val- 
ues reported [19]. However the slightly higher values of 
nic = 1.5 GeV and mi, = 4.7 GeV are also used in the lit- 
erature. Changes in quark mass will affect the fragmen- 
tation functions. In the scheme of our calculation, the 
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TABLE II: Total cross section in pb for triply heavy baryons in possible c and b quark fragmentation for various scales at the 
LHC with ^/s = 14TeV where the kinematical cuts of pt > 10 GeV and|j/| < 1 are imposed. Note that the cross sections are 
calculated in two groups of scales, {fi = Hr/2, 2/iH and 4/x_r) for lighter c -^ 0.ccc, c -^ Qccb and b —> riccb and (/i — hr, Z^r 
and Q^r) for heavier c -^ ilcbb, b -^ ilcbb and 6 —> Qbbb states. The ratios a(Q -^ Q,)/a{Q) are given in the last column. They 
are calculated at /^ = 2/i_R. 



Process of 








Cross Section 


[pb] 






The Ratio 


Production 


Mii/2 


IJ.R 


2/iii 


S^iR 




4:fiR 


6fJ.R 


a{Q ^ n)/a{Q) 


C — > Qccc 


301.88 




306.99 






307.59 




1.20 X 10"' 


C -^ Qccb 


26.58 




30.03 






29.88 




1.18 X 10"* 


b -^ Qccb 


2153.08 




2155.31 






1723.80 




3.61 X lO"'^ 


b -^ ritbb 




6.34 


6.38 


5.77 






8.40 


1.36 X 10"** 


b -^ i^cbb 




50.30 


34.77 


47.78 






52.34 


7.43 X 10-** 


c -^ Qm 




1.14 


1.38 


1.47 






1.49 


5.41 X 10"^" 



fragmentation functions inversely depend on quark mass 
squared. Therefore increase in quark mass will decrease 
the probabilities. The other quantity which may depend 
on quark mass is the baryon decay constant. However the 
later is not much clear in the case of triply heavy baryons. 
Taking the explicit mass dependence of our fragmenta- 
tion functions, we have obtained 18 percent decrease in 
the cross sections in average, when we use the above men- 
tioned higher values. 

There is no data on the baryon decay constant. Theo- 



retically one may solve the Schrodinger like equation to 
obtain the wave function at the origin for these compos- 
ite particles with heavy constituents and then relate the 
wave function at the origin to the baryon decay constant. 
We have avoided this procedure because of theoretical 
uncertainties instead have chosen fs = 0.25 GeV on phe- 
nomenological grounds. The final quantity of interest is 
the color factor. We have calculated this quantity us- 
ing the simple color line counting rule and have obtained 
Cf = 7/6 for our propose. 
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